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t— I . ABSTRACT 

I"*"""- ' We discuss the environmental dependence of galaxy evolution based on deep panoramic 

imaging of two distant clusters, RXJ0152. 7-1257 at z = 0.83 and CL0016+1609 at 
z = 0.55, taken with Suprime-Cam on the Subaru Telescope as part of the PISCES 
project. By combining with the SDSS data as a local counterpart for comparison, we 
construct a large sample of galaxies that spans wide ranges in environment, time, and 
• stellar mass (or luminosity) . This allows us to conduct systematic and statistical anal- 

yses of the photometric properties of galaxies based on the colour-density diagrams, 
q . colour-magnitude relations, and luminosity functions. We find that colours of galax- 

' ies, especially those of faint galaxies (My > M v + 1), change from blue to red at a 

\ break density as we go to denser regions. This trend is observed at all redshifts in our 

Ci • sample. Based on local and global densities of galaxies, we classify three environments: 

field, groups, and clusters, and look into the environmental dependence of galaxies in 
. i-H ! detail. In particular, we quantify how the colour-magnitude relation is built-up as a 

■ function of environment. We show that the bright-end of the cluster colour-magnitude 

^ ' relation is already built at z — 0.83, while the faint-end is possibly still in the process 

of build-up. In contrast to this, the bright-end of the field colour-magnitude relation 
has been vigorously built all the way down to the present-day and the build-up at the 
faint-end has not started yet. A possible interpretation of these results is that galaxies 
evolve in the 'down-sizing' fashion. That is, massive galaxies complete their star for- 
mation first and the truncation of star formation is propagated to smaller objects as 
time progresses. This trend is likely to depend on environment since the build-up of 
the colour-magnitude relation is delayed in lower-density environments. Therefore, we 
may suggest that the evolution of galaxies took place earliest in massive galaxies and 
in high density regions, and it is delayed in less massive galaxies and in lower density 
regions. Further studies are, however, obviously needed to confirm the observed trends 
and establish the 'down-sizing' picture. 

Key words: galaxies: fundamental parameters — galaxies: evolution — galaxies: 
luminosity function, mass function — galaxies: clusters: individual CL0016+1609 — 
galaxies: clusters: individual RXJ0152.7— 1357 
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1 INTRODUCTION 

Intensive studies on galaxy properties, such as star forma- 
tion rates and morphology, have significantly improved our 
understanding of galaxies in the Universe. It is, however, 
still unclear how galaxies evolve over the Hubble time. This 



is due to the complex nature of galaxy properties; galaxy 
properties depend not only on time, but also on environ- 
ment and mass (luminosity). These three axes are related to 
one another and they characterize galaxy evolution. Despite 
the obvious importance, however, galaxy properties along 
these three axes have not been viewed simultaneously and 
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systematically, due to the limited previous data sets none of 
which can cover all of these three axes. Based on the wide- 
field Subaru data and the large SDSS data, we present in 
this paper a comprehensive study of star formation activity 
of galaxies along the three axes. 

It is now well known that galaxy properties depend on 
environment in which galaxies reside. This environmental 
dependence of gala x y pro perties was first quantitatively 
studied by IPressleJ l)l980l) , who showed the morphology- 
density relation based on 55 nearby clusters. There is 
a clear trend that early-type galaxies are preferentially 
located in high density regions, while late-type galaxies 
tend to be located in lower density regions. Following 
this work, intensive studies on environmental depen- 
dence of galaxy properties have been carried out, and 
strengthen ed or ex tende d the Dressler's result (e.g., 
Postman fc Gelled Il984t IWhitmore Gilmore. fe Jonea 
199 1 IBaloeh et al.lll997t Dressier et alJll997t iBaloeh et all 
19981: lHashimoto et alJ Il99l iLubin et alJ Il998t 
Qjte^Post^nan^ ^^^^to|" | l99St IPostman. Lubin. fc Okel 



199sl: Ivan Dokkum etaf ll99St IBaloeh et all Il999t 



Poggianti et al.1 Il999l: ILubin et all l2000l ICouch et al 
2001 : iKodama et alJ |2001al: IPostman. Lubin fc Oke 
200ll: iLewis et alJ 1200a ILubin Oke fc Postman 1200 " 
Blanton et alj|2003al: iGomez et alJ | 200l IGoto et al 112003. 



Hogg et alJ 120031 iTreu et alJ 12003: IBaloeh et alJ l2004allbl 
Smith et al.lbo04t ITanaka et al]l2004l) . 

Galaxy properties are also known to depend on 
mass (luminosity) of gal axies. Recent l y con ducted large 
surveys, such as 2dF dColless et alJ 12001) and SDSS 
iYork et alJ I2000TI . revealed that gal axy properties show 



strong bimodality in t heir distribution ( Stoateva^^lJ|20ryil 
Blanton et a i] |2003al : iKauffmann et alJ l2003t IBaloeh et alJ 



2004al:lKauffmann et alJl2004l : ITanaka et al.ll2004) . That 



there are two distinct populations: red early-type galax- 
ies and blue late-type galaxies. This bimodality is found 
to be a strong function of mass of galaxies in the sense 
that massive galaxies tend to be red early-type galaxies, 
while less massive galaxies tend to be blue late - type galaxies 
12001 



ijK auffmann et al 
200 



iBalogh et alJ 120044 iBaldrv et all 
iBell et alJl2004t IKauffmann et alJl2004Tl 

Galaxies evolve with time, of course, as a nat- 
ural consequ ence of stellar e v olution . It is now 20 
years since iButcher fc Oemlerl lll984h presented their 
startling result that the fraction of blue galaxies 
in clusters increases with look-back time or red- 
shift. This B-O effect seems to be confirmed by 



later studies (ICouch & Sharoles 


Il987l iRakos & SchombertJ 


1991 lEllingson et al.l 1200 ll 


Kodama & Bowerl l2001bl 


Margoniner et al. 


20011 iNakata et alJ 1200 ll iFairlev et all 


2001 IGoto et al 


2003bl iTran et alJ 12004 but see also 


Fairlev et alJl2002l 


Andreon. Lobo. & Iovindl2004l). This ef- 



feet suggests that star formation activities in cluster galaxies 
change with time, in the sense that the fraction of star form- 
ing galaxies decreases with time. 

These three axes, that is, environment, mass, and 
time, must be closely related to one another. For exam- 
ple, the morphology-density relation is found to evolve 
llDressler et alJ Il997t iFasano et all 12001 ITreu et al.1 l200l 
Smit h et al .120041 but see also lAndreonll998h . ITa naka et alJ 
( 2004) showed that environmental dependence of star forma- 
tion activity and morphology of galaxies change as a func- 



tion of galaxy luminosity (mass). To untie this complexity, 
we base our analyses on panoramic imaging data of two high 
redshift clusters, CL0015.9+16 at z = 0.55 and RXJ0152.7- 
13 at z = 0.83, and the SDSS data. The former data are ob- 
tained as part of our on-going Subaru distant cluster project 
called PISCES. Our data span wide ranges in environment, 
time, and stellar mass, and give us an unique opportunity to 
investigate star formation activity along all the three axes 
simultaneously for the first time. 

The structure of this paper is as follows. In fJ21 we briefly 
review our project PISCES and our data. Also we describe 
the data of the local Universe that we use. Then we move 
on to the procedure we adopt to eliminate foreground and 
background contamination in J|] Before presenting main re- 
sults, we summarize various photometric and environmental 
quantities used in this paper in We examine environ- 
mental dependence of galaxy star formation in Based 
on results obtained in f0 we re-define environments in 
Colour-magnitude diagrams and luminosity functions in var- 
ious environments are shown in Sj7|and|H| respectively. In SjSJ 
we discuss the implications of our results on the formation 
and evolution of galaxies. Finally, a summary is given in H10I 
Most of the important conclusions in this paper are drawn 
from the results presented in §5.1, §7, and §9. Readers in- 
terested only in our primary conclusions can go to these 
sections directly after brief reading of §2-4. 

Throughout this paper, we assume a flat Universe of 
Q M = 0.3, tt A = 0.7 and H = 70kms~ 1 Mpc" 1 . We use the 
AB magnitude system for observed quantities and the Vega- 
referred system for rest-frame ones. We use the following 
abbreviation: CMD for colour-magnitude diagram, CMR for 
colour-magnitude relation, and LF for luminosity function. 



2 DATA 

2.1 PISCES Project 

We are conducting a systematic study of cluster evolution 
based on panoramic multi-band imaging with Suprime-Cam 
and optical multi-slit spectroscopy with FOCAS on Sub- 
aru. In this section, we briefly review the project Panoramic 
Imaging and Spectroscopy of Cluster Evolution with Subaru 
(PISCES). The reader should refer to lKodama et al.l (12001) 
for further details. 

The primary aims of this project are two- fold: (1) to 
map out large scale structures around distant clusters to 
trace the cluster assembly history and (2) to look into galaxy 
properties in detail as a function of environment along the 
structures in order to directly identify the environmental 
effects acting on galaxies during their assembly to higher 
density regions. The unique feature of this project is its 
wide- field coverage by tak ing advantage of the Suprime- 
Cam (iMivazaki et all 2002). which provides a 34' x 27' field 
of view. Therefore, the PISCES will provide an opportu- 
nity to link an evolutionary path between galaxies in the 
local Universe and high redshift counterparts over a wide 
range in environment. Also, the depth of PISCES, reach- 
ing down to M* + 4 at z ~ 1 with an 8-m telescope, 
will shed light on the nature of faint galaxies at high red- 
shifts, which probably show quite different properties com- 
pared with massive galaxies, as seen in the local Universe 
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(e.g., iBaldrv et all 12004 iKauffmann et alJl2004h . As part 
of this on-going project, we observed CL0015. 9+1609 and 
RXJ0152.7-1357 in September 2003. 

The cluster CL0015.9+1609 (CL0016 for short) is 
one of the most extensively studied galaxy clusters, and 
it has be en observed in various wav elengt h ranges: X- 
ray (e.g., IWorrall fc Birkinshawl l2003h . UV ferown et alJ 
|2QQ_fJ), optical (e.g., iDjessLer^ ^dllTQQflh . and submm (e.g., 
IZemcov et ail booaTlkocI dl98 lb suggested the existence 
of a large-scale structure around the cluster based on the 
photographic photometry, and later, the cluster was found 
to ha ve companion clu sters (iHughes. Birkinshaw. fc Huchral 
Il995l: IConnollv et alJll995 lHugh.es fc Birkinshawlll998li. A 
clear CMR is seen llEllis et alJ Il997t iDahlen et al.ll2004) . 
and the cluster is known to have a very low blue fraction 
jButcher fc Oemlerlll98l . 

The cluster RXJ0152. 7-1357 (RXJ0153 for short) is 
one of the most X-ray luminous distant (z > 0.7) clus- 
ters known. The cluster was discovere d independently in 
the ROSAT Deep Cluster Survey (RDCS: lRo"sati et alll99gl) 
and in the Wide Angle RO SAT Point ed Survey (WARPS; 
IScharf et aHll997t Ebeling et al] l2000ft . Later, it was also 
detected in the Serendipitous Hiq h-redshift Arc hival ROSAT 
Cluster (SHARC) survey jRomer et aljl2000h . Since these 
discoveries, the cluster has been the subj ect of BeppoSAX, 
Chandra and XMM-Newton observations ( Delia, Ceca et al] 
I2OOOI: iMauehan et~ai1 l2003i: IJones et al.l ho03|. O bserva- 
tions of the Sunyaev- Zeldovich effect <Jov et al.l200l. spcc- 
troscopic follow-up jDemarco et alJ l200^ ~lHo^^iCT|et_al 
2004; Ijctrgensen et all2004f ). near-IR imagi ng fellfs^^oneiL 
200#). and the weak lensing anal ysis llHuo^tHnTT2"o"o4i 



2UU4), and t ne we ak lensmg analysis HJluo et alJ Vz 
Umetsu et all 120051: IJee et al] [2004b have also been 



per- 
formed. All these studies show that RXJ0153 is indeed a 
very massive cluster at z = 0.83: the bolometric X-ray lumi- 
nosity of > 1 x 10 45 ergs s _1 , and the total dynamical mass 
of ~ 1 x 10 15 M p). Based on the ACS weak lensing analysis, 
IJee et al] (j2004T) suggested that the previous cluster mass 
estimates may be an over-estimation. 

Below, we briefly summarize our observation of these 
two clusters. De tails of the obs e rvatio n and data reduction 
are described in lKodama et all <l2005b . 

The observing conditions were excellent during the 
nights with a typical seeing size of 0".6. CL0016 was 
observed in BVRi'z' and RXJ0153 in VRi'z'. Exposure 
times and limiting magnitudes are shown in Table 
Object detection is pe rformed using SExtractor (v. 2. 3. 2; 
Bertin fc Arnoutslll996l) . Objects in CL0016 are i'-band se- 
lected and those in RXJ0153 are z'-band selected. It should 
be noted that the i'-band for CL0016 and the z'-band for 
RXJ0153 correspond to almost the same rest-frame wave- 
length range, and therefore we have little difference in selec- 
tion effects between the two clusters. We use MAG_AUTO 
for total magnitudes and 2" diameter aperture magnitudes 
for colours. Magnitude zero-points determined by the photo- 
metric standard stars taken during the same observing run 
are fou nd to be offset (< 0.1 mag.) with respect to stellar 
SEDs of lGunn fc Strvkerl lll983b . We shift the zero-points so 
as to match with Gunn & Stryker stars. Star-galaxy separa- 
tion is performed on the basis of FWHM vs. total magnitude 
diagrams. 

In what follows, we restrict ourselves to i' < 24.5 galax- 
ies in CL0016 and z < 25.0 galaxies in RXJ0153. These 



Cluster 


Filter 


Exposure Time (min.) 


Limiting Magnitude 


CL0016 


B 


90 


26.9 




V 


96 


26.2 




R 


64 


26.0 




i> 


60 


25.9 




z' 


47.5 


24.6 


RXJ0153 


V 


120 


26.7 




H 


116 


26.5 




i' 


75 


26.1 




z! 


77 


25.0 



Table 1. A list of exposure times and limiting magnitudes (AB 
system). Limiting magnitudes are shown as a 5<r limit measured 
in a 2" aperture. 



magnitude cuts ensure that we are not affected by incom- 
pleteness effects. 



2.2 SDSS 

As a local counterpart of the two high-z samples, we use the 
data from the Sloan Digital Sky Survey (SDSS:lYork et al] 
l200rl IStoughton et alJ 120021: lAbazaiian et alJ l200l 12004b 



i 



The SDSS observes one quarter of the sky both photometri- 
cally and spectroscopically. The imaging sur vey is performed 
in five optical bands, u, q, r, i, and z ([Fukugita et aF 
| l99dlGunn et alJll99a iHogg et aUl200ltlSmith et alJbo'o " 
iPier et al]l2003t) . The spectroscopic survey is done with a 
pair of double fiber-fed spectrographs which covers 3800A - 
9200A. Each fiber subtends 3" on the sky. Since a fiber can- 
not be placed closer than 55" to a nearby fiber due to me- 
chanical constraints, the tiling algorithm has been developed 
to red uce the number of unobserved objects jBlanton et alJ 
l2003bb . The overall completeness of the spectroscopic survey 
is expected to be over 90%. 

We use the the public d ata of the second data release 
(DR2: lAto : zaj iaii_et : ^Jj2jW^ . Galaxies in the Main Galaxy 
Sample Tsfaauss et aljl2002b are used here. The sample cov- 
ers 2627 square degrees of the sky. We extract galaxies at 
0.005 < z < 0.065. We perform a volume correction to a 
flux limited sample taking into account large-scale struc- 
tures, instead of making a volume-limited sample, so that 
we can statistically reach as deep as the high-z samples 
in terms of magnitude relative to the characteristic mag- 
nitude. Details of the volume correction is described in Ap- 
pendix |A"| We use Pe trosian magnitudes iPetrosianl Il976l : 
IStoughton et all2002Tl for total magnitudes and model mag- 
nitudes for col ours. All these quantities are Galacti c extinc- 
tion corrected iSchlegel. Finkbeiner. fc Davidllliih . and k- 
corrected using the code of Blanton et al. (2003c, v3_2). 

Main galaxies in the SDSS are r-band selected 
jStrauss et alJl2002h . while those in RXJ0153 and CL0016 
are selected in the rest-frame ~ (?-band. To minimize selec- 
tion effects, we apply the magnitude cut of g < 18 to our 
SDSS sample so that the sample mimics a (/-selected one at 
the cost of reducing the number of sample galaxies. This cut 
leaves 41695 galaxies. In summary, the galaxies in RXJ0153, 
CL0016 and SDSS are all selected at similar wavelengths in 
the rest-frame nearly corresponding to the g-band. 
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3 CONTAMINATION SUBTRACTION 

Since our two clusters from PISCES lie at high redshifts 
and our imaging is deep, galaxies at the cluster redshifts are 
heavily contaminated by foreground and background galax- 
ies. Thus, in order to study galaxies at the cluster redshifts, 
it is essential to eliminate the contamination. Our strat- 
egy for the contamination subtraction is two-fold. Firstly, 
we apply photometric redshift tech nique to largely elimi- 
nate fore- /backgro und galaxies fe.g.. iKodama et alJ feoOla: 
iNakata et al.l20(W I. Although photometric redshift is a pow- 
erful tool, the contamination remains at a non-negligible 
level. We therefore statistically subtr act the remaining con- 
tamination on the b asis of CMDs ijKodama et alJ 2001a; 
iPimbblet et al]l2002ft . Each procedure is described in the 
following subsections. 



3.1 Photometric Redshifts 

We apply the photome tric redshift code of 
IKodama. Bell, fc Bowerl il999l) to the photometric cata- 
logue of the PISCES clusters (IKodama et al]l2005l) . The 
p hotometric redshift ut il izes th e population synthesis model 
of IKodama fc Arimotol {j.997). Star formation histories of 
model templates are described by a combination of the 
elliptical galaxy model with a very short time scale (O.lGyr) 
of star formation and the disk model with a much longer 
time scale of r = 5 Gyr (star formation rate oc e~ t ^ T ). The 
models are construc ted so as to reproduce the o bserved 
colours of galaxies llKodama. Bell, fc Bowerl 11999"). It is 
found that the observed colours of the CMR are slightly 
offset compared with the model colours. We shifted the 
model zero-points systematically so that the model colours 
of the CMR match with the observed colours at the cluster 
redshifts. These shifts are required to calibrate and improve 
our photometric redshifts. 

To assess the accuracy of our photometric redshifts, we 
compare the photometric redshifts with spectroscopically 
determined redshifts as shown in Figure ^ For RX J0153, 
we use the spectroscopic sample of iJorgensen et alJ (J2004). 
Our photometric redshifts are fairly good for cluster galax- 
ies. But, at z < 0.5, the accuracy is rather poor since we 
lack U and _B-band data. Excluding the most deviant galaxy 
at z apec = 0.745, the mean and median of z v hoto — z apec 
arc +0.0067 and +0.0093, respectively. The standard devia- 
tions around the mean and median are found to be 0.045 and 
0.047. For CL0016, we compile spectro s copic redshift data o f 
Hughes. Birkinshaw. fc Huchral <1995h. iMunn et all |l997tl . 
Hughes fc Birkinshawl <1998ft . and iDressler et al.l ll99Sli . 
Excluding the two most deviant galaxies at z apec ~ 0.4, 
the mean and median of z p hoto ™ z apec are +0.0036 and 
—0.0118, respectively. The standard deviations around the 
mean and median are 0.055 and 0.057. The figures demon- 
strate that our photometric redshifts are fairly accurate. 
We note, however, that the spectroscopic samples are het- 
erogeneous, and the face values quoted above should not 
be over-interpreted. Note as well that the photometric er- 
rors of red galaxies at our magnitude limit are typically 
a(V) = 0.1, a(R) = 0.06, cr(i') = 0.04, cr(z') = 0.06 
in RXJ0153 and o(B) = 0.09, a(V) = 0.05, a{R) = 
0.02, cr(i') = 0.02, a(z') = 0.04 in CL0016 (errors of blue 
galaxies are smaller than these values). Thus, photo-z should 



work with a reasonabl y good accuracy (|Az| < 0.1) even at 
our magnitude limits iKodama. Bell, fc Bowerlll999l) . 

Using cluster galaxies only (0.78 < z ap( , c < 0.88 for 
RXJ0153 and 0.50 < z spec < 0.60 for CL0016), we exam- 
ine the accuracy of photometric redshifts as a function of 
galaxy colour. For RXJ0153, we cannot identify any bias in 
the photometric redshifts, although we lack galaxies having 
intermediate colours (R— z ~ 1.4). We find, however, that 
we tend to underestimate redshifts for galaxies with interme- 
diate colours (V-i' ~ 1.5) by Az ~ -0.1 in CL0016. A pos- 
sible reason for this colour dependence is that slightly bluer 
galaxies at cluster redshift than those on the CMR tend to be 
regarded as red galaxies at slightl y lower redshifts because 
of th e colour-redshift degeneracy iKodama. Bell, fc BoweJ 
1999). But, we need more spectroscopic data to fully ad- 
dress this issue. Our photometry covers a similar rest-frame 
wavelength regime for both RXJ0153 and CL0016, and we 
expect that we also tend to underestimate redshifts of galax- 
ies with intermediate colours in RXJ0153. In both samples, 
photometric redshifts of red galaxies are found to be fairly 
accurate, though a small offset z p \ loto — z apec — +0.02 is seen. 

Galaxies outside of a certain photometric redshift range 
around the cluster redshift are regarded as fore-/background 
galaxies. The choice of the redshift range is a trade off be- 
tween the completeness of cluster galaxies and the contam- 
ination of fore-/background galaxies. If we adopt a small 
redshift range, the contamination will be reduced, but a 
price to be paid is a selection bias towards red galaxies. 
If a wide redshift range is adopted, the selection bias will 
be reduced at the cost of increasing an amount of the con- 
tamination. Since we are interested in galaxy properties, es- 
pecially colours, we prefer to construct an unbiased sample. 
For this purpose, we adopt z c i — 0.12 < z p hot < z ci + 0.05, 
where z c i is a cluster redshift. This selection criteria will 
eliminate the colour selection bias well, while maintaining 
an amount of the contamination to be minimal (Figure 0. 
To be specific, we adopt 0.42 < z phot < 0.60 for CL0016 
and 0.71 < z p hot < 0.88 for RXJ0153. Such an asymmetric 
redshift ranges are taken because of the possible asymmetric 
photo-z error distributions (see Fig. . 

It should be noted that these redshift r a nges are dif- 
ferent from those adopted in IKodama et alJ ll2005f) . They 
adopted narrower redshift ranges to enhance large-scale 
structures by tracing, primarily, red galaxies. 

3.2 Statistical Contamination Subtraction 

Even though the photometric redshift is effective to reduce 
the fore-/background contamination, contamination still re- 
mains at a non-negligible level within our photometric red- 
shift ranges. We construct a control field sample and statis- 
tically subtract such remaining contamination on the basis 
of the galaxy distribution on CMDs. Note that we do not 
apply any contamination subtraction in the SDSS sample 
since galaxies are spectroscopically observed. 

Here we describe only the essence of the procedure of 
the statistical subtraction. Sec Appendix IB1 for details. We 
adopt a modified procedure of IKodama fc Bowerl {2001b) 
and lPimbblet et alJ i2002fl . In brief, the distribution of con- 
trol field galaxies on the CMD is used as a probability map 
of the contamination. The control field sample is defined as 
low-density regions in each field of RXJ0153 and CL0016, as 
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Figure 1. Left: Photometric redshifts plotted against spectrally determined redshifts for RXJ0153. In the inset, we show differences 
between photometric redshifts and spectroscopic redshifts as a function of R — z' colour, which is close to the rest-frame U — V colour, 
using galaxies at 0.78 < z spec < 0.88. Right: Same as the left panel but for CL0016. In the inset, we use V — i' colour and galaxies at 
0.5 < z spec < 0.6. 



shown later. The choice of the control field is arbitrary, but 
our results do not strongly rely on a particular choice of the 
control field. We select two separated regions in each field 
to minimize cosmic variance. The average galaxy density in 
the selected regions is denoted as E co „t ro ;. The subtraction 
of the remaining contamination is done on the CMD using 
a Monte-Carlo method. We statistically subtract the galaxy 
distribution on the CMD of the control field from that of 
the target field. 

3.3 Concerns about the Contamination 
Subtraction 

Although we carefully subtract the contamination, uncer- 
tainties arising from the subtraction cannot be ignored. In 
this subsection, we briefly address the robustness of our con- 
clusions presented below. 

In the following sections, we discuss the fraction of red 
galaxies. As shown above, photometric redshifts of red galax- 
ies are more accurate than those of blue galaxies, and we 
may tend to miss blue galaxies. Therefore, the real fraction 
of red galaxies would be smaller than we observe. This will 
strengthen our conclusions, the observed density dependence 
of galaxy colours and evolutionary trends. We cannot, how- 
ever, evaluate the amount of missing blue galaxies with the 
data at hand. 

Errors in the statistical subtraction is also a concern 
(e.g., over/under subtraction). However, taking advantage of 
the wide field of view, we can take wide areas for our control 
field sample, namely ~ 270 arcmin 2 and ~ 210 arcmin 2 in 
RXJ0153 and CL0016, respectively. Therefore cosmic vari- 
ance must be averaged to some extent and is not expected to 
be a serious problem. We repeated Monte-Carlo runs of sta- 
tistical subtraction many times and confirmed that trends 
we discuss later are seen in most realization. In other words, 



we discuss secure results only. Although the robustness of 
our contamination subtraction must be confirmed later spec- 
troscopically, it is unlikely that our results significantly suffer 
from the uncertainties in the contamination subtraction. 



4 DEFINITIONS OF ENVIRONMENTAL 
PARAMETERS AND DERIVATION OF 
REST-FRAME QUANTITIES AND STELLAR 
MASS 

Before presenting our results, we summarize various photo- 
metric and environmental parameters used in this paper. 
We begin by introducing environmental parameters used 
to characterize local and global environments. We then de- 
scribe galaxy properties such as magnitude, colour, and stel- 
lar mass. 



4.1 Environmental Parameters 

4-1.1 Local Density 

In this paper, we use nearest-neighbor density to charac- 
terize environment mainly because it is a frequently used 
indicator and comparisons with other studies can be made 
directly. In RXJ0153 and CL0016, all the galaxies in the 
selected photometric redshift range are projected onto the 
redshift of the main cluster, and density is estimated from 
the distance to the 10-th nearest galaxy from the galaxy of 
interest and is denoted as Eioth- We use a circular aperture 
in the density calculation. This density is called local density 
hereafter. It should be noted that local density is actually 
a surface density. Galaxies that reach the edge of our field 
of view before finding their 10-th nearest galaxies are not 
used in the analysis since the local density of such galaxies 
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is not correctly estimated. Local density for the two high- z 
clusters is calculated in both physical and comoving scales. 

In the SPSS, local density is estimated in a similar 
manner as iBaloeh et alJ <l2004bf) . In brief, galaxies within 
±1000 km s -1 in the line-of-sight velocity space from the 
galaxy of interest are projected onto the redshift of the cen- 
tral galaxy, and local density is defined from the distance to 
5-th nearest neighbor. When counting galaxies, we use only 
those brighter than My = —19.5, which is volume-limited 
out to z = 0.065 (our maximum redshift range). This den- 
sity is denoted as Esth- Although we use the distance to the 
5-th nearest neighbor, our results are essentially unchanged 
if we use the 10-th nearest neighbor as used for RXJ0153 
and CL0016. Since the contamination of fore-/background 
galaxies is very small in the SDSS sample, we prefer to adopt 
a smaller number so that density represents more 'local' en- 
vironments. The fiber-collision problem may affect our den- 
sity estimates in high-d ensity environment s. However, the 
effect is not significant dTanaka et al.l l20041. Local density 
for the SDSS is calculated in the physical scale only. There is, 
however, little difference between the physical and comoving 
density since galaxies lie at very low redshift s. 

4-1-2 Global Density 

Global density is defined as the surface galaxy density in 
a fixed radius of 2 Mpc around the central galaxy. Galax- 
ies are projected onto the cluster redshift (RXJ0153 and 
CL0016) or onto the redshift of the galaxy (SDSS) in ques- 
tion, in just the same manner as local density. Since we aim 
to characterize global environment, global density is eval- 
uated as a comoving density. By combining local density 
with global density, we can separate poor groups from rich 
clusters quantitatively. The effectiveness of this method is 
described in Appendix[0] In what follows, global density is 
denoted as E gtobat . 

4.2 Magnitudes, Colours, and Stellar Masses 

4-2.1 Rest-frame Magnitudes and Colours 

For RXJ0153 and CL0016, we derive rest-frame absolute 
U-band magnitude and U — V colour from the observed 
magnitudes and colours us ing the model templates of th e 
photometric redshift code (iKodama. Bell, fc Bowerlll999l) . 
The conversion applied in the observed quantities to derive 
the rest-frame quantities is small. Within the selected red- 
shift ranges, the effect of galaxy evolution is estimated to 
be AM V < 0.2 and A(U - V) < 0.05 for both RXJ0153 
and CL0016. As for the SDSS sample, we estimate the rest- 
frame V and U — V using the code by Blanton et al. (2003c, 
v3_2). We recall that, for conventional reasons, we use the 
Vega-referred system in the rest-frame V magnitude and the 
U — V colour. Table [5] summarizes the limiting magnitudes 
for the three samples. 

4-2.2 Stellar Masses 

We derive approximate stellar masses of galaxies. For each 
sample, the stellar mass to light ratio (M*/L) and M* are 
derived from the model templates used in the photometric 
redshift code. Our L is defined at the rest-frame ~ U-band, 



sample 


Mv,lim 




SDSS 


-17.5 


4 x 10 9 


CL0016 


-18.0 


5 x 10 9 


RXJ0153 


-18.0 


4 x 10 9 



Table 2. A list of the V limiting magnitudes (Vega) and the 
limiting stellar masses of the three samples. 



and so the estimates of M* are affected by on-going/near- 
past star formation activities. It is expected that our stellar 
mass estimates are relatively accurate for red galaxies, since 
red galaxies are not actively forming stars. However, stellar 
masses of blue galaxies cannot be reliably determined. We 
find that M* / L ratios span a factor of ~ 4 depending on the 
colour of galaxies, and we can get only rough estimates in 
stellar mass although we correct for the mass-to-light ratio 
based on the colours (SED fitting). Moreover, stellar mass 
is a model-dependent quantity, since mass-to-light ratio de- 
pends on stellar init ial mass function (IMF). O ur M* esti- 
mate is based on the IKodama fc Arimotd (1199711 population 
synthesis model, and we assume the IMF of x = 1.10 for the 
elliptical models and x = 1.35 for the disk models with the 
mass range of O.IMq — 6OA/0. The limiting stellar masses 
that we can trace completely are shown in Table [5] 



5 COLOUR-DENSITY RELATIONS 

We apply photometric redshift technique and discover large- 
scale structures around bot h RXJ0153 and CL00 16 clus- 
ters. Details are described in IKodama et all i2005l) . Due to 
the wide-field coverage of the Suprime-Cam, we obtain a 
wide variety of environments, i.e., sparse fields, poor groups, 
and rich clusters. Here we examine the relationship between 
galaxy colours (star formation rates) and environments. We 
stick to local environments in this section. Effects of global 
environments are examined later. First, we characterize en- 
vironment by surface galaxy density (local density). Next, 
environment is defined on the basis of surface mass density. 

5.1 Dependence on Surface Galaxy Density 

Based on d a ta fro m lar ge surveys of the l ocal Universe, 
iLewis et all (l2002f) and iGomez et all J2003I) showed that 
galaxy star formation begins to decline sharply, in a sta- 
tistical sense, at a certain local density. In what follows, this 
sharp decline is referred to as 'break', and the density where 
'break' is seen is referred to as break density. Non-star- 
forming galaxies dominate regions above the break density, 
whereas star-forming galaxies are the dominant population 
below the break den sity. Based on the data from the SDSS, 
iTanaka et alJ (120041) showed that the break is seen only for 
galaxies fainter than M* + 1, and brighter galaxies show no 
clear break in the local Universe. Following this work, we ex- 
amine the environmental dependence of star formation for 
bright and faint galaxies separately. Galaxies brighter than 
M v + 1 are defined as bright galaxies, and those fainter than 
that limit are defined as faint galaxies (the value of My at 
each redshift is given in §S|. 

Figure [5] shows our results. In the figure, the U — V 



The Build-up of the Colour-Magnitude Relation 7 



colour is corrected for the slope of the CMR so that galaxies 
on the CMR have the same colour as that of a My galaxy 
and is denoted as (U — V) CO rr (i.e., CMR is transformed into 
a horizontal sequence). The faint galaxies (solid line) show 
a break, prominent change in their (U — V) corr colours, at 
the densities marked by the dot-dashed lines in the figures. 
It is interesting that the break in galaxy colours is seen at 
all redshifts examined here. In contrast, the bright galaxies 
(dashed line) do not show such a strong break, especially 
the median lines. The (U — V) COT r colours of bright galaxies 
are systematically redder than those of faint galaxies at any 
local density. 

In order to quantify the break feature, we measure the 
slope of the colour-density relation using the median lines 
in Fig. |21 as a function of local density for bright and faint 
galaxies separately. The result is shown in Fig. 01 The faint 
galaxies show a strong change in slope at the densities shown 
by the dot-dashed lines, while bright galaxies show a much 
weaker change there. Here we define the break density at 
which we see the strong colour change in faint galaxies, as 
marked in Figs. [SI and |3] 

We change the threshold magnitude used for separat- 
ing bright/faint galaxies and check how the presence of the 
break density of each cluster depends on the magnitude of 
galaxies. It is found that galaxies brighter than My + 1 do 
not show a prominent break, while those fainter than that 
limit show a break at the same density. 

There are two possible effects that cause the break. 
One is that the fraction of red galaxies relative to blue 
galaxies begins to increase above the break density. The 
other is that blue galaxies become systematically redder 
(but bluer than galaxies on the CMR) above the break den- 
sity. These possibilities are investigated in the right panels 
of each plot in Figure [5] Red and blue galaxies are sepa- 
rated at (U - V)cmr - 0.15. In RXJ0153 and CL0016, the 
(U — V) CO rr of blue galaxies does not change with density, 
while the fraction of red galaxies strongly changes with den- 
sity. This means that the break is caused by the change 
in the population fraction of red galaxies. As for the SDSS 
plot, the (U — V) C orr of faint blue galaxies becomes system- 
atically redder above the break density. On the other hand, 
the (U — V) C orr of bright blue galaxies does not strongly 
change with density. The fraction of red galaxies is found to 
strongly depend on density. Therefore, we conclude that, in 
the SDSS plot, the break is caused by both effects. 

The break densities are used to define environments in 
JS] We note that a small change in the break density has no 
significant effect on our results. Note as well that the break 
density at z = 0.55 and z = 083 is 3-5 times higher than 
the control field density, and thus the uncertainty in the 
statistical contamination subtraction is not a concern. As 
described in §3.3, photometric redshifts may miss a fraction 
of blue galaxies. This will strengthen the observed break 
since the break is primarily driven by the change in the 
population fraction of red galaxies relative to blue galaxies. 
Further discussion on the break density of each cluster is 
made in 89.2. 



5.2 Dependence on Surface Mass Density 

iGrav et al.l (120041) first presented the environmental varia- 
tion of galaxy colours as a function of surface mass density. 
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Figure 3. The slope of the colour-density relation [ A(t7 — 
V)/ A log S ] as a function of local density along the median loci 
in Fig. |21 The top/bottom ticks are comoving/physical densities 
as in Fig.H The panels show RXJ0153, CL0016, and SDSS (from 
left to right). The filled/open symbols show bright/faint galax- 
ies, respectively. The vertical dot-dashed line means the break 
density. The errors are estimated from the bootstrap resampling. 



Inspired by this work, we discuss colours of galaxies as a 
function of surface mass density in this subsection. This in- 
vestigation is particularly interesting since surface galaxy 
density, on which our analysis in the previous subsection 
is based, represents density of luminous matter around the 
cluster, while the surface mass density estimated via the 
weak-lensing analysis represents total mass including dark 
matter. These two densities do not necessarily agree with 
each other, and comparisons between the dependence on 
galaxy density and that on mass density will give us a hint 
of physical mechanisms that affect galaxy properties. 

The weak lensi ng mass reconstructi on of RXJ0153 is 
described in detail in lUmetsu et alJ i2005ll . The lensing con- 
vergence k, is related to the surface mass density E K by 
k(9) = E re (0)/E^. ff crit . As a fiducial value, lUmetsu et alJ 
J2005D adopted E** rit = 3.1 x lO 15 /i 7 oM Mp C - 2 . 

We plot in the left panel of Figure 0] the surface galaxy 
density against the lensing convergence n(0). A positive cor- 
relation is found between k and galaxy surface density, es- 
pecially at high densities with /t £ 0.1. On the other hand, 
no clear correlation can be seen in the low density regime. 

In the right panel, the (U — V) cor r colour is plotted 
against k. Although it is not as clear as seen in Figure |5J 
there is a hint of a break in the (U — V) corr colour at k ~ 0.1. 
This threshold c orresponds t o E K ~ 3 x 10 14 M Q /i 70 Mpc" 2 
in physical units. IGrav et alJ ||2004) found a similar break at 
the surface mass density of E K ~ 3.6 x 10 14 M ftroMpc" 2 , 
which is consistent with our estimate. However, this appar- 
ent threshold density of k ~ 0.1 is comparable to the rms 
noise level in the reconstructed ft map, a K ~ 0.10. There- 
fore, the underlying mass density threshold can be smaller 
than what we obtained, k ~ 0.1. It is therefore premature to 
say if galaxy properties are more strongly related to galaxy 
density than to mass density. We note that lJee et alJ (l2004ll 
recently reported that distribution of galaxies, mass, and in- 
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Figure 2. The colour-density relation in RXJ0153 (top-left), CL0016 (top-right) , and SDSS (bottom). These plots show one realization 
of our Monte-Carlo run for the statistical contamination subtraction. Note that no contamination subtraction is performed in the SDSS 
plot. Left panel in each plot: The (U — V) C orr colour plotted against local density. The lines show the median and 25-th percentile of 
the distribution of bright (dashed line) and faint (solid line) galaxies as noted in the RXJ0153 panel. The median line is associated with 
bootstrap 90% intervals as shown by the error bars. Densities are expressed in the comoving density (top ticks) and in the physical 
density (bottom ticks) and are contamination corrected (i.e., shifted leftward by S contro j). The vertical lines show the break density of 
each cluster. The arrow indicates ^control, where a half of galaxies are statistically subtracted. In the SDSS plot, local density is shown 
in a physical scale only. For clarity, one tenth of all the SDSS galaxies are randomly selected and plotted. Each bin contains 100 bright 
/ 200 faint galaxies in RXJ0153 and CL0016, and 1000 bright / 2000 faint galaxies in the SDSS plot. Top-right panel in each plot: The 
(U — V) C orr colour plotted against local density for blue [U — V < (U — V)cmr — 0.15] galaxies. The lines show the median of the 
distribution of bright/faint galaxies as noted in the RXJ0153 panels. The associated error bars are the bootstrap 90% intervals. The 
horizontal line means (U — V)cmr — 0.15. The vertical line shows the break density of each cluster. Bottom-right panel in each plot: 
The fraction of red galaxies plotted against local density. The meanings of the lines are given in the RXJ0153 panels. The errors are 
based on the Poisson statistics. 



tracluster medium are all different in this cluster suggesting 
on-going cluster merger. 



Although we cannot draw a firm conclusion on this anal- 
ysis, this is potentially an interesting way of investigating 
environmental dependence of galaxy properties. If, for ex- 
ample, galaxy-galaxy interactions are the main driver of the 
environmental dependence, we expect to see stronger de- 
pendence on galaxy density than on mass density. On the 
other hand, if interactions with cluster potential is the main 
driver, we expect a stronger dependence on mass density. 



6 DEFINITIONS OF FIELD, GROUP, AND 
CLUSTER ENVIRONMENTS 

In order to address galaxy properties as functions of envi- 
ronment and time, the definition of the environment must be 
quantitative and applicable at all redshifts. For this reason, 
we define environments by galaxy density. For analyses in 
the following sections, we define three environments: field, 
groups, and clusters as shown in Figures |S] and0 

We do not examine galaxies in lower density regions 
than ^control because most of them are considered to be 
fore-/background galaxies. We find in the previous section 
that red galaxies dominate environments denser than the 
break density. Motivated by this, we use the break density 
of each cluster to define environments. Galaxies in higher 
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Figure 4. Left: Relationship between the surface galaxy den- 
sity (local density) and normalized surface mass density (k) in 
RXJ0153. The lines show the median and quartiles (25% and 
75%) of the distribution. The horizontal dashed line means the 
break surface galaxy density. The arrow indicates k = 0.1 where 
S/N of the mass map is unity. Right: The (U — V) C orr colour 
plotted against k. The meanings of the lines are the same as Fig 
|21 The arrow indicates mass density of S/N = 1. 



density regions than E contro ;, but in lower density regions 
than the break density (Ebreak) are considered to be in the 
field environment (i.e., E co , ltro ( < Y>i oca i < Inbreak)- 

Galaxies in higher density regions than E(, rea fc belong 
to groups and clusters. We now aim to investigate the 
effects of global environment: are there any differences 
between group galaxies and cluster galaxies? This question 
is particularly interesting since different physical mecha- 
nisms are effective in different enviro nments. For example, 
ram-pressure stripping of cold g as llGu nn fc Gottl |l972t 
[Abadi. Moore, fc Bowei] Il999t fOuilis. Moore, fc Bowed 
1200(f) and harassment iMoore et alJ Il996l. ll99Sh are ex- 
pected to play a role only in r ich clusters. Low- velocity 
galaxy-galaxy interactions (e.g., iMihos fc Hernq uist 1996) 
are expected to be effective in galaxy groups. Strangulation 
(which is often referred to as suffocation, starvation or 
halo gas stripping) is considered to be effective both in 
groups and clusters llLa rson^Tins^^ 

iDiaferio et alJ 



Balo gh. Navarro, fc Morris! 1 2000 ; 



Okamoto fc Nagashimal l2003 ) . 



1980; 



2001; 



Global density defined in H4.1.2I is found to work well 
in separating groups from clusters. Global density traces 
galaxy density over a large scale, and thus global density 
is a good measure of richness of galaxy systems. We define 
globally denser environment than the break density as clus- 
ters and globally less dense environment as groups. That is, 
clusters are defined by E ioc(l i > E breafc and S s i ob(l! > Y> brea k, 
and groups by E; oca! > E breafc and E siobai < T, break . As 
illustrated in Figures |S| and [7] our separation of envi- 
ronments is reasonable. We show in each plot the virial ra- 
dius of the cluster. The v irial radius of RXJ0153 is taken 
from lMaughan et alJ (I2003T) . and those of CL0016 and SDSS 
are evaluated from velocity dispe rsions using the recipe of 
ICarlberg. Yee. fc Ellingsonl Jl997|) . As seen in the plots, the 
break density corresponds to the outskirts of clusters. This 
is quantitatively consistent with that seen in the local Uni- 



Environment 


Definition 






Field 


^control 


< ^local < ^brcak 




Group 


^local > 


Sbrcak 


and S g i b a i 


Sbrcak 


Cluster 


^local > 


^brcak 


and S global 


> ^break 



Table 3. Definitions of environments. S C ontroli ^locab ^brenk, 
and Sgi bai are control field density, local density (S 5th or Ej^^), 
break density, and global density. 



verse where the break density corresponds to one to two 
R vir llGomez et alJl2003l : iTanaka et alj|2004t) . We note that 
the break density also corresponds to a density typical of 
isolated groups. 

We are aware that, since we define environment based 
on galaxy colours, we may obtain 'biased' environmental de- 
pendencies of galaxy colours. For example, since we define 
cluster environment where red galaxies are abundant, clus- 
ters are, by definition, dominated by red galaxies. However, 
the fact that the break density corresponds to the outskirts 
of clusters is a good justification of our definition. There- 
fore, we consider that results presented below are not biased 
products of the environment definition. 

To sum up this section, we tabulate the definitions of 
environments in Tabled] Based on these environments, we 
look into CMDs and LFs of galaxies in the following sections. 



7 COLOUR-MAGNITUDE DIAGRAMS 

Gala xies in clusters are know n to show a tight CMR (e.g., 

This relation is observed 



up to z = 


1 and even beyond llKodama et alJ 


1998; 


Stanford, Eiscnhardt, & Dickinson 1998; Nakata et alJ 


2001: Blakcslcc et alJ 12004 iDe Lucia et all 


2004; 


Lidman et alJ 


2004|) ■ In this section, we aim to 


inves- 



tigate the CMDs with particular attention to the group and 
field environments, which have not been intensively studied 
yet especially at high redshifts. 

Figure |H| shows the rest-frame CMDs. We estimate the 
slope and scatter of the CMRs based on an iterative 2a- 
clipping least squares fit and the results are shown in Fig- 
ures l9land !10l Note that the CMR slopes are measured from 
galaxies brighter than My- +2 which show tight relationship, 
while the scatter around the CMR is measured for galaxies 
down to our magnitude limits. Errors are estimated from 
the bootstrap resampling of the input catalogs, and the sta- 
tistical field subtraction is performed in each run. Since our 
photometric redshift cuts have some ranges (Az ~ 0.18), a 
projection effect is a concern which may apparently enhance 
the intrinsic colour scatter. However, the colour difference in 
passively evolving galaxies within these redshift ranges are 
only 0.05 magnitude at most for both CL0016 and RXJ0153 
llKodama et alJl998l) , therefore this effect is small compared 
to the amount of scatter we discuss here (>0.1 mag.). We 
see five very interesting trends in these figures. 

Firstly, we find probable onset of the build-up of the 
CMR in field regions. In the field environment in RXJ0153 
at z — 0.83, we cannot identify any clear CMR, although a 
clump of bright red galaxies can be seen. In fact, the scatter 
around the CMR is very large (Fig.|5Jl, being consistent with 
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Figure 5. Galaxy distribution in RXJ0153. North is up. The filled rectangles and open circles show cluster and group galaxies. The 
large dots are field galaxies. Galaxies in the lower density regions than S contro ; are shown by the small dots. Galaxies in the large dashed 
rectangles are used as control field galaxies, and they are used in the statistical contamination subtraction procedure. Galaxies too close 
to the edge of our field of view are not plotted for they have n o local density estimates. The top and right ticks show the comoving scales 
in unit of Mpc. The circle shows R v i r (Mauehan ct al. 200^). No statistical contamination subtraction is applied in this plot. 



no CMR. Interestingly, however, a clear relation is seen in 
the same environment in CL0016 at z = 0.55 and also in 
SDSS at z = particularly at the bright-end. We find that 
the scatter around the CMR at the bright-end decreases 
from z = 0.83 down to z = 0, while the scatter at the faint- 
end does not decrease clearly. That is, the bright-end of the 
CMR is built-up with time, while there is no clear CMR at 
the faint-end even at z = 0. 

Secondly, we find on-going build-up of the CMR in 
groups. Groups in RXJ0153 show a CMR, but we find that 
the CMR is getting weak at Mv > —20. There are a num- 
ber of red galaxies at Mv > —20, but they do not form a 
tight relation. The scatter around the CMR shown in Fig. 
[5] reflects this visual impression. At lower redshifts, group 
galaxies show a clear relation down to the magnitude limit. 
This leads us to suggest that we are witnessing the build- 
up of the CMR in groups. An emerging picture is that the 
CMR grows from the bright-end to the faint-end, and not 
the opposite way. 

Thirdly, no such on-going build-up of the CMR is seen 
in clusters. We cannot see dramatic growth of the CMR 
in clusters. Cluster galaxies show a clear CMR down to the 
magnitude limit at all redshifts considered here and the scat- 
ter around the CMR is already small at high redshifts. 

One may suspect that the observed build-up of the 
CMR is an artifact: e.g., photometric redshifts miss a frac- 



tion of red galaxies in RXJ0153, and it mimics the CMR 
build-up. But, this is not the case. The accuracy of pho- 
tometric redshifts does not depend on environment. Ac- 
cordingly, if there were a CMR in the field environments 
of RXJ0I53, it would have been found because we see the 
clear CMR in the cluster environments. The same argument 
can be applied to the faint-end of the CMR in groups. It 
should be noted that our group environment is a composite 
of several individual groups, and hence a variation in group 
properties is not a major concern. Note as well that photo- 
metric redshifts are reliable for red galaxies (see Figure 0. 
Therefore, we suggest that the observed build-up is real. 

Fourthly, the slope of the CMR does not change with 
environment. The slope of the CMR in the field environ- 
ment is found to be similar to those in groups and clusters. 
This is quantified in Figure 1101 There is no convincing ev- 
idence that the CMR slope depends on environment. Note 
that galaxy colours are measured within different aperture 
sizes at different redshifts: Petrosian aperture in SDSS and 
2" aperture in CL00 16 and RXJ0153. CMR slop es change 
with aperture sizes jBower. Lucev. fc Ellisl IT992j) . and we 
do not discuss the evolution of the CMR slope here. Al- 
though CMR seems to be built-up at different epochs in dif- 
ferent environments, this similarity of the CMR slope may 
be expected if the slope is primarily caused by the m etal- 
licity effect rather than the age effect iKodama &: Arimotol 
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or equivalently, per unit magnitude, 



Il997l: IStanford. Eisenhardt. fc Dickinson! Il998lh The ob- 
served CMR build-up suggests that the slope of CMR does 
not change since the formation epoch of the CMR. 

Fifthly, we confirm bimodality in galaxy colours. 
Galaxy properties are known to have strong bimodal 
distribution in the l ocal Universe (|S tra±eva_el_alJ |200jJ; 
Blanton et"al] l2003al: iKauffmann et al] 120031: iBaldrv et all 
20041: iBalogh et al l l2004bl IKauffmann et al] l2004t 
Tanaka et alJl2004h . This bimodality is seen up to z ~ 1 
]Bell e^u*Tl2C)04l) . We confirm the bimodality although 
it becomes less clearly seen at higher redshifts in all the 
environments. The bimodality is particularly noticeable 
in the field regions of CL0016 and SDSS. High density 
environments lack blue galaxies. 

The highlight of this CMR analysis is that we observe 
the build-up of the CMR. It seems that an evolutionary 
stage of the CMR build-up is different in different environ- 
ments: the cluster CMR is built first, and the CMRs of the 
group and field regions are built later on. Another interest- 
ing implication is that the bright-end of the CMR appears 
first, and the faint-end is filled up later on. We will further 
quantify and discuss the observed build-up in §9.3. 



8 LUMINOSITY FUNCTIONS 

Luminosity function (LF) is one of the most fundamental 
measures of galaxy properties. This section studies LFs as 
functions of environmen t and time . A LF is fitted by the 
Schechter function iSchechterlll976f) of the form 

^ L )dL=<j>*(^y B^L^di^j, (i) 



4>{M)dM = 0.4 In 10 x ^iq-o^m-m^^+V) 

exp [_i -°- 4 ( M - M *)] d M. (2) 

There are three parameters in the Schechter function: 
the normalization factor <f>*, the characteristic luminos- 
ity/magnitude L* JM* , and the faint-end slope a. Best-fit 
parameters are searched via the conventional x2-minimizing 
statistics. Galaxies in RXJ0153 and CL0016 are binned into 
0.5 magnitude steps, and those in the SDSS are into 0.25 — 
0.5 magnitude steps. First, we examine total (=red+blue) 
LFs at each redshift. Then we look into red/blue LFs in 
different environments and at different redshifts. 



8.1 Total Luminosity Functions 

Based on the photo-z selected samples of RXJ0153 and 
CL0016 and the spec-z sample of SDSS, the total LF of 
galaxies at each redshift is constructed and shown in Fig- 
ure Note that no statistical contamination subtraction 
is performed here. 

The Schechter function gives a good fit (xt ~ 1) to the 
total LF of for RXJ0153 and CL0016, but it is not a good fit 
for SDSS (xl ~ 11). The observed SDSS LF deviates from 
the Schechter function at the faint-end, and this deviation 
decreases the overall goodness of fit. The deviation is possi- 
bly due to increasing contribution of dwarf galaxies. We fit 
the Schechter function using galaxies with Mv < —18.5 and 
obtain M£ = -21.12 and a = -1.01 with xl ~ 4. We re- 
call that the My derived in Fig. 1111 were used in separating 
bright/faint galaxies in |5] Note that a small error in M v 
has little effect on the results obtained in that section. 

Figure 1111 clearly shows that galaxies fade with time: 
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Figure 8. The rest-frame CMDs (U - V vs. M v ) in RXJ0153 (left), CL0016 (middle), and SDSS (right) are plotted. In each plot, 
the panels show different environments, namely, field, group, and cluster. Small dots mean statistically subtracted galaxies (except the 
SDS S plot). The solid line represent the CMR shifted blueward by 0.15 magnitude. The CMR is based on the model prediction of 
iKodama fe Arimotol ll997l) . This solid line separate red and blue galaxies used in the following analysis. The vertical dashed line is the 
magnitude limit, and the slanted dashed line shows the 5<r limiting colour. The arrows show My + 1, which is used to separate bright 
galaxies from faint galaxies in JF] In the cluster plot, the curved dashed lines show loci of stellar masses of 10 10 and 10 11 Mq. Note 
that, in the SDSS plot, the volume correction is applied and the corrected galaxies are artificially plotted. Note as well that no statistical 
contamination subtraction is performed in the SDSS plot. 



My = -22.05 at z = 0.83, My = -21.61 at z = 0.55, and 
M v = -21.24 at z = 0. This observed fading, AM V = 0.81 
mag. from z — 0.83 to z = and AMy = 0.37 mag. from z = 
0.55 to z = 0, is consistent with a passive evolution model 
jKodama fc Arimotdll997t: z f = 5), AM V = 0.80 mag. and 
0.55 mag., respectively, within the errors. This suggests that 
My is primarily determined by passively evolving galaxies. 
In our cosmology, z = 0.83 and z = 0.55 correspond to a 
look-back time of 7.0 Gyr and 5.4 Gyr, respectively. 

The faint-end slope a also seems to evolve. The number 
of faint galaxies relative to bright galaxies increases with 
time, a = -0.94 at z = 0.83 and a = -1.12 at z = 0. But 
the local value should be taken with caution since a is found 
to depend on the magnitude range involved in the Schechter 
fit. Contribution of dwarf galaxies is likely to be significant 
at the faint-end. We therefore do not try to draw any firm 
conclusion on the evolution of a. 

In RXJ0153 and CL0016 fields, we have rich clusters 
and a number of groups, and the contribution of clusters and 
groups to the total LF is larger compared with the SDSS LF. 
In the SDSS sample, the field, group, and cluster galaxies 
comprise 80%, 13%, and 7% of the total number of galax- 
ies. We scale the relative fraction of field, group, and cluster 
galaxies in RXJ0153 and CL0016 to the SDSS values, and 
find that the derived Schechter parameters show only a small 
change: AMy = 0.1 and Act = 0.02. Thus, a different en- 



vironmental mix of galaxies does not strongly change our 
results. 



8.2 Luminosity Functions of Red/Blue Galaxies 

The fact that galaxy properties have strong bimodality in 
their distribution motivates us to investigate LFs of red and 
blue galaxies separately. We define red galaxies as those hav- 
ing U — V > (U — V)cmr — 0.15. Blue galaxies are defined as 
those bluer than this limit. This definition is illustrated in 
Figure [S] It should be noted that, in the following, the sta- 
tistical contamination subtraction is performed on the LF 
bin - LF bin basis. This is different from the Monte-Carlo 
approach that we adopt in the previous sections. Our re- 
sults are presented in Figure [T2*l A general trend is that red 
galaxies have a decreasing or flat faint-end slope (though 
exceptions can be found), while blue galaxies have an in- 
creasing slope at any redshift. Some SDSS data points at 
the faint-end deviate from the best-fit Schechter function. 
This is because the SDSS data at the faint-end have small 
weights in the Schechter fitting due to the statistical nature 
of the volume correction. 

Now, we focus on the LFs of red galaxies. At z = 0, 
the faint-end slope becomes g radually less steep in denser 
environment. iHogg et all (120031) reported, based on SDSS 
data, that faint red galaxies preferentially populate in high 
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Figure 11. Total LFs in RXJ0153 (top), CL0016 (middle), 
and SDSS (bottom). Note that no contamination subtraction is 
performed here. The left panels show the error ellipse for each 
LF. The solid and dashed contours show xS = xt best ^ ant ^ 
= xl ibest + 2, where xl, best is the xl of the best fit. The right 
panels show LFs along with the best-fit Schechter functions. The 
error bars are based on the Poisson statistics only. 



densi ty environments. A similar conclusion was reached by 
iDe Prooris et alj l)2003ft based on the 2dF data. It is likely 
that red faint galaxies are selectively located in high density 
environments. A similar trend can be seen in RXJ0153, but 
it is not statisti c ally s ig nificant given t he lar g e error ellipses. 
iKaiisawa et alJ fcOOOft. iNakata et all fcOOlft. IrCodama et alJ 



J2004D . iToft et all i2004t >. and IDe Lucia "et alT 12004) re- 
ported the deficit of red faint galaxies in high redshift 
(z > 0.7) clusters compared with local clusters. Although 
we cannot confirm the trend in the Schechter parameter, we 
do see the the deficit of red faint galaxies in the giant-to- 
dwarf ratio, which we will discuss in the next section. As 
for blue LFs, there is no convincing evidence that blue LFs 
strongly depend on environment and time. 



9 DISCUSSION 

9.1 Possible Peculiarity of CL0016 

Galaxies in CL0016 show somewhat exceptional properties 
compared with other samples. For example, as shown in Fig. 
1121 the field environment in CL0016 has a strange red LF, 
and the group environment has an increasing faint-end slope 
for red galaxies. One may find that these LFs look different 
compared with other LFs. The fraction of red galaxies is the 
highest in all environments as we see later. 

There is a concern that this peculiarity may arise from 
errors in the statistical subtraction of contamination. We 
therefore check how the results for CL0016 can change by 
changing our choice of field regions used in the subtraction. 
We take either of the two field regions indicated by rectan- 
gles in Fig. [5] separately, rather than combining these two 
fields. We then compare the resulting CMDs and LFs for 
CL0016 after the statistical subtraction. This effectively cor- 
responds to changing the field density by factor of ~ 2. Yet, 



we find no evidence that the results change strongly. There- 
fore, we suggest that the peculiarity of the CL0016 cluster 
is intrinsic. CL0016 is probably one of the oldest systems in 
the Uni verse where red faint galax ies are already abundant. 
Indeed, iButcher fc Oemlerl (^984) illustrated the unusually 
low blue fraction of this cluster for its redshift. Our environ- 
ments are defined on the basis of local and global (2Mpc) 
density. It might be the case that even larger-scale envi- 
ronments (e.g . , ~ 10 Mpc) can aff ect the galaxy prope rties 
jBaloeh et al.ll2004al: but see also iBlanton et alJl2004l) . In 
the following, we regard CL0016 as an exceptional sample. 



9.2 Break Density 

The break density corresponds to the outskirts of 
galaxy clusters (one to two R V i r ) and isolated groups. 
This means that groups and clusters are dominated by 
red galaxies and supports the idea that cluster-specific 
mechanisms have not played major roles i n transform- 
ing ga lax y properties jKodama et aljfeoOlal; i Baloeh et alJ 
l2004allbt foe Propris et alj |2004 iTanaka et al] 12004 but 
see also iFuiital l2004t iFuiita. fc Gotol 12004ft . Rather, 
mechanisms such as low-velocity ga laxy-galaxy interac- 
tions (e.g.. iMih os fc Herncmisd | l996ft and strangulation 
jLarson. Tinslev. fc Caldwelllll980ft" remain as strong candi- 
dates of the driver of the environmental dependence. 

The fact that faint galaxies (My > My + 1) show the 
clear break suggests that physical mechanisms actually work 
on faint galaxies in high-density environments. Then, why 
do not bright galaxies (My < My + 1) show any prominent 
break? Of course, bright galaxies are generally red and a 
break in their colours, if any, would be difficult to see. But, 
at least in the SDSS plot in Figure^ this is not the case. No 
break is seen even in the 25-th percentile of bright galaxies. 
This implies that the evolutionary path of bright galaxies is 
different from that of faint galaxies in such a way that the 
evolution of faint galaxies is strongly related to groups and 
clusters, while that of bright galaxies is not strongly related 
jTanaka et al.ll2004ft . 

Environmental dependence of galaxy properties 
is determined by a priori effects (initial conditions) 
and a posteriori effects (environmental effects) . Re- 
cent near-IR studies such as K20 Survey iCimatti et alJ 



Daddi et alJ 120021: iFontanal 



|2004ft. 



FIRES 



( Franxet^ljj2p03; van Dokk um et al.l l2003: Ru dnick et alJ 
120031: ISchreiber et all l2004ft ancl part of GOODS 
jGiwalisc^^^^J]200?lD^dre t al J 120041 : iMoustakas et all 
l20o"4T ~ TsOTnBrvmeet^rtJ 12004ft" revealed the existence of 



massive red galaxies at 2 > 1. Although a significant 
fraction of such red galaxies are dusty starbursts (e.g., 
iMivazaki et alJ 12003ft . massive non-star- forming galaxies, 
which are presumably precursors of present-day ellipticals, 
do exist. Since their colours match with passive evolution, 
their properties are expected to be largely determined by 
a priori effects and subsequent environmental effects are 
not very important. That is, their star formation rate is 
already low before environmental mechanisms play a role. 
These massive galaxies would evolve to bright galaxies in 
our definition, and would explain the trends observed in 
On the other hand, the deficit of red faint galaxies is 
a function of environment, in a way that red faint galax- 
ies preferentially populate in high-density environments. 
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Figure 12. Rest-frame LFs in various environments in RXJ0153 (top-left), CL0016 (top-right), and the SDSS (bottom). The left panels 
in each plot show the error ellipses (%J = \^ best + 1) of the Schechter fits. The solid and dashed lines represent red and blue galaxies 
defined in Figure|S] The filled points and open squares show the best-fit parameters for red and blue LFs, respectively. The right panels 
show LFs. The filled and open symbols represent red and blue galaxies. The solid and dashed lines respectively show best-fit Schechter 
functions for red and blue galaxies. The statistical contamination subtraction is carried out here. 



Therefore, we suggest that properties of bright galaxies are 
largely determined by a priori effects, while those of faint 
galaxies are largely determined by a posteriori effects. 

9.3 The Build-up of the CMR 

In <J7| we observed the build-up of the CMR. To further 
quantify the build-up, we base our analysis on stellar masses 
of galaxies. If galaxies stop their star formation, they will be 
fainter (AMy ~ 1) while keeping their stellar masses nearly 
unchanged. Since our stellar mass estimates are not very ac- 
curate, we cannot examine a detailed shape of a stellar mass 
function. Instead, we investigate the giant-to-dwarf number 
ratio (g/d). Giants and dwarves are defined as those hav- 
ing log 10 (M,/M@) > 10.6 and 9.7 < log lo (Af,/M ) < 10.6, 
respectively. Note that log 10 (M./M©) = 10.6 corresponds 
to My = —20.5, —20.3, and —19.8 for red galaxies at 
z = 0.83, 0.55, and 0, respectively. 

Since our estimates of stellar masses are primarily based 
on the the rest-frame V-band magnitudes, they are af- 
fected by recent /on- going star formation activities. In fact, 
as shown by the iso-mass curve in Fig. |HJ the l/-band mag- 



nitude can change by 1.5 mag. for the same stellar mass 
between passively evolving galaxies and the constantly star 
forming galaxies. This , which is then translated to the vari- 
ation in mass-to-light ratio by a factor of ~ 4. This can be 
viewed as a solid upper limit of the uncertainties in M* in 
a relative sense, since we correct for such variation in mass- 
to-light ratio by applying SED fitting when deriving stellar 
masses (note that absolute stellar masses depend on other 
factors such as the stellar initial mass function). Moreover, 
since we mainly discuss the red galaxies, the actual varia- 
tion in mass-to-light ratio must be much smaller (less than 
a factor of 2). In what follows, we focus on the red galaxies 
only. We separately discuss field, group, and cluster environ- 
ments. 

Field Environment: Figure tTSI shows g/d in three dif- 
ferent environments and at three different redshifts. The pa- 
rameter g/d clearly depends on both environment and time. 
A general trend is that g/d is the largest in the field environ- 
ment at any redshifts, suggesting that red faint galaxies is 
relatively rare in the field. This is consistent with our finding 
in Fig. |S] that the faint-end of the field CMR is not clear. 
The field g/d ratio is the largest in SDSS. This is likely to 
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be driven by the build-up of the bright-end of the field CMR 
(i.e., the fraction of giants increases). 

Group Environment: The g/d ratio in groups is the 
largest in RXJ0153 and it decreases at low redshifts. This 
decrease in g/d may reflect the build-up of the group CMR 
at the faint-end: the fraction of faint red galaxies increases. 

Cluster Environment: The cluster g/d also shows the 
decrease at low redshifts. Although we could not see promi- 
nent build-up of the cluster CMR in 33 the g/d evolution 
may suggest the faint-end of the cluster CMR is still under 
construction even at z — 0.83. 

It seems that the bright-end of the CMR appears first, 
and the faint-end is filled up later on. A likely scenario for 
this build-up is that blue galaxies stop their star formation 
and fade to settle down onto the CMR, and the trunca- 
tion of star formation starts from bright (massive) galaxies. 
This scenario involves suppression of star formation activ- 
ity in blue galaxies. How do blue galaxies in groups and 
field environments fade and settle exactly on the CMR as 
that of cluster galaxies? We recall that we found no ev- 
idence for environmental dependence of the slope of the 
CM R. If the CMR is a product of the mass-metallicity rela- 
tion iKodama fc Arimotdll99'ft) . then blue galaxies have to 
'know' the metallicities of red galaxies they settle onto. Since 
blue galaxies in the field and group environments should 
follow quite different star formation histories from those of 
red galaxies in clusters, one may expect that their metal- 
licities are quite different. This is not the case, however. 
An important point here is that we are considering rela- 
tively massive (log 10 M* > 9.7) galaxi es. The cosmic star 
formation rate declines at z < 1 (e .g., |Mjdjm_j^t_ajJ 1996; 
iMadau. Pozzetti. fc DickinsorJl99alFuiita et all2003ft . and 
therefore their major episode of star formation took place 
at higher redshifts. Subsequent star formation does not 
strongly enrich metals in galaxies, and metals locked in their 
stars do not grow s ignificantly after the major star forma- 
tion jTmslevlll98Ch . Therefore, once galaxies mature, the 
epoch when they stop their star formation is not a major 
concern from a chemical point of view. Whenever they stop 
their star formation, they will settle onto the CMR. 

To further quantify the build-up of the CMR, we plot 
in Figure ITU the fraction of red galaxies as functions of en- 
vironment, stellar mass, and time. As described in §3.3, we 
cannot deny the possibility that we miss a fraction of blue 
galaxies in RXJ0153 and CL0016. Therefore, the real red 
fractions in RXJ0153 and CL0016 would be smaller than 
those shown in the plot. We discuss the red fraction in the 
field, group, and cluster environments separately as follows. 

Field Environment: The massive galaxies show an in- 
crease in the red fraction, which is consistent with the build- 
up of the bright-end of the field CMR. On the other hand, 
the least massive galaxies show a decrease. However, this 
may be because we tend to miss a fraction of blue galaxies 
in RXJ0153 and CL0016 due to our photo-z selection. 

Group Environment: Due to the relatively large er- 
rors, the red fraction of group galaxies is consistent with 
being unchanged with redshift. 

Cluster Environment: Within the errors, red fraction 
is consistent with being almost constant over the time under 
study. This is in agreement with the trend that no significant 
CMR build-up is seen in clusters within the redshift range 
explored. There is an hint of evolution, however, at the least 



massive bin, where red fraction may increase from z = 0.83 
to lower redshifts. 

The overall trend is that the red fraction is the lowest in 
the field environment and higher in groups and clusters. But, 
in any environments and at any redshifts, the red fraction is 
higher for more massive galaxies. That is, the red fraction is 
higher in denser environment and for more massive galaxies. 

9.4 Implications for Galaxy Evolution 

From Figs.HSlandnn we consider that an evolutionary stage 
of the CMR build-up is different in different environments. 
The bright-end of the cluster CMR is already built by z = 
0.83, and only the faint-end shows a possible evolution at 
z < 0.83. On the contrary, in the field regions, the bright- 
end of the CMR is being vigorously built, while the build-up 
at the faint-end has not started yet. 

A possible interpretation of these results is that strong 
evoluti on has occurred sin ce z ~ 1 in a 'down-sizing' way. 
It was lCowie et alJ dl99rj) who first pointed out the down- 
sizing galaxy evolution. At high redshifts, massive galaxies 
actively form stars. At lower redshifts, massive galaxies show 
less active star formation and the main population of active 
star formation moves to less massive galaxies. The bright- 
end of the CMR is formed first as massive galaxies stop their 
star formation, and the build-up proceeds to the faint-end 
as less massive galaxies stop their star formation. Our re- 
sults suggest that the evolutionary stage of this down-sizing 
depends on environment. The evolution of massive cluster 
galaxies is almost completed by z = 0.83. In the field envi- 
ronment, the evolution of massive galaxies is strong, while 
that of less massive galaxies is found to be weak. There- 
fore, it seems that the main population that shows strong 
evolution is shifted to higher mass galaxies in lower density 
environments. 

One of the drivers of the environmental dependence of 
the down-sizing could be initial conditions of galaxy forma- 
tion (a priori effects). Galaxies are formed earlier in higher 
density peaks of the initial density fluctuation of the Uni- 
verse. Thus, galaxies in clusters are naturally at an advanced 
stage of galaxy evolution compared with those in lower- 
density environments. This may explain the environmental 
dependence of the build-up of the CMR. We consider how- 
ever that the environmental dependence of the down-sizing 
effect is not solely caused by the a priori effects, but envi- 
ronmental (or a posteriori) effects should contribute signif- 
icantly. Effects that suppress star formation activities are 
strong in high-density environments (3SJ, and they acceler- 
ate the build-up of the CMR. 

We observed the build-up of the CMR at the faint-end 
in group environments. This should mean either that faint 
galaxies are still actively forming stars, or that less mas- 
sive galaxies are not fully formed yet. We cannot, however, 
discriminate these possibilities with the data in hand. To 
do this, we need to estimate stellar masses of faint blue 
galaxies accurately. Further discussion on the driver of the 
CMR build-up awaits more accurate stellar mass estimates 
by near-infrared data. 

Detailed observations of early-type galaxies suggest that 
the typical luminosity-weighted age of fiel d early- type galax- 
ies is younger than cluster galaxies (e.g., |Kuntschner et alJ 
120021: ICebhardt et al.l2003f) . lThomas et al.l(l2004l) examined 



The Build-up of the Colour-Magnitude Relation 17 



1 0.5 



redshift 
1 0.5 1 



0.5 



$ RXJ0153 
+ CL0016 
4 SDSS 



ID 

d 




field 



group 



cluster 



Figure 13. The giant-to-dwarf ratio of red galaxies plotted 
against environments. The meanings of the lines are shown in 
the figure. The errors are based on the Poisson statistics. 



nearby early-type galaxies and suggested that formation of 
early-type galaxies is the earliest in high density environ- 
ments and delayed by ~ 2 Gyr in low density environments, 
and formation of massive galaxies pred ates that of less mas- 
sive galaxies. Utiliz ing near-IR data, iFeulnerl j2004) and 
lJuneau et alJ l|2004f) reported that the epoch of major star 
formation took place at higher redshifts for more massive 
galaxies, and star formation is more extended in time for 
less massive galaxies. All these studies lend support to the 
down-sizing picture. 

We close our discussion by noting some caveats on our 
results. Our results are based only on two high-z clusters, 
which may not be a typical cluster at each redshift (par- 
ticularly CL0016). To reach a firm conclusion, we need to 
observe more clusters at various redshifts. Also, errors in the 
photometric redshifts and stellar mass estimates remain as 
a concern. Near-infrared (e.g., if-band) data are required to 
improve them. Spectroscopic data are clearly needed to fur- 
ther address the effectiveness of the photometric redshifts 
and assess errors in the statistical contamination subtrac- 
tion. We hope to overcome these uncertainties in our future 
paper. 



10 SUMMARY AND CONCLUSIONS 

We began this paper by introducing the three axes on which 
galaxy properties strongly depend, namely, environment, 
stellar mass, and time. We found that the star formation 
activity in galaxies is indeed dependent on all of these three 
axes and galaxies follow complicated evolutionary paths. 

We conducted multi-band imaging of two high-z clus- 
ters, RXJ0153 at z = 0.83 and CL0016 at z = 0.55, with 
Suprime-Cam on Subaru. These Subaru data were combined 
with the SDSS data (z = 0), and we carried out statistical 
analyses of galaxy properties as functions of environment, 
mass, and time. We examined the colour- density relations, 
colour- magnitude diagrams (CMDs), and luminosity func- 
tions (LFs) of galaxies. 

First, we applied our photometric redshift tech- 
nique to RXJ0153 and CL0016 fields to largely eliminate 
fore-/background contamination and discovered large-scale 
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Figure 14. Fraction of red galaxies plotted against the look-back 
time (the corresponding redshift is shown in the top ticks) for the 
three stellar mass bins. The meanings of the lines are shown in 
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structures surrounding the ma in clusters. Details are de- 
scribed in lKodama et all i2005h . 

Then we examined the relationship between galaxy 
colours and environments. It was found that galaxy colours 
abruptly change at the break density at any redshifts con- 
sidered in this paper. Faint galaxies (My > M v + 1) show a 
prominent break, while bright galaxies (My < M v + 1) do 
not show such a strong break. Based on the break density, 
we defined three environments: field, group, and cluster. 

We examined CMDs in the three environments. The 
highlight of the CMD analysis was to show how the CMR is 
built up as functions of time, environment and mass. There 
is no clear CMR in the field regions in RXJ0153, while a clear 
CMR can be seen in the same environment in CL0016 and 
SDSS, particularly at the bright-end. In groups of RXJ0153, 
a clear CMR is visible only at the bright-end (My < —20). 
The relation is found to extend down to our magnitude limit 
in groups of CL0016 and SDSS. Clusters have clear CMRs 
down to our magnitude limits at all redshifts considered 
here, with a possible remaining activity at the faint-end. 
These trends are quantified by the scatter of the CMRs. 

The build-up of the CMR was quantified by the giant- 
to-dwarf ratio (g/d) and the red fraction. From z = 0.83 
to z = 0, the field g/d increases and the group g/d de- 
creases, suggesting the CMR build-up at the bright-end in 
the field and at the faint-end in the group. The cluster g/d 
also decreases. This may imply that the faint-end of the clus- 
ter CMR was still under construction at high redshifts. We 
found that the red fraction is higher in denser environment 
and for more massive galaxies. 

As a possible interpretation of these results, we sug- 
gested that galaxies evolve in the down-sizing way. That 
is, the main populations that host active star formation is 
shifted from massive galaxies to less massive galaxies with 
time. It is likely that an evolutionary stage of the down- 
sizing depends on environment. Cluster galaxies evolve most 
rapidly, and group and field galaxies follow. All in all, it 
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seems that galaxy evolution proceeds from massive galaxies 
to less massive galaxies, and from high density environments 
to low density environments. In order to confirm this trend, 
however, further studies are obviously needed. We need to 
increase the cluster sample and also to perform follow-up ob- 
servations such as NIR imaging and spectroscopic surveys. 
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APPENDIX A: DETAILS OF THE VOLUME 
CORRECTION 

In a flux-limited sample, intrinsically faint galaxies are ob- 
served only at low redshifts, while intrinsically bright galax- 
ies can be observed at higher redshifts. For example, in 
our SDSS sample, we can observe galaxies brighter than 
Mv, vega < —19.5 at all redshifts, while My, V ega < — 18 can 
be observed only at z < 0.033. We therefore need to correct 
for this incompleteness when we discuss galaxy properties. 
A simple way to do this is to give heavier statistical weights 
to intrinsically fainter galaxies according to their redshifts. 
A flux-limited sample can mimic a volume-limited one in 
this way. A statistical weight can be computed as 

W^i — Vsurvey / Vi,max , (-^-1) 

where w% is a statistical weight of i-th galaxy, V 3U rvey is the 
volume contained within our redshift range (0.005 < z < 
0.065), and Vi t max is the maximum volume over which the 
i-th galaxy could be observed (e.g., in case of Mv = —18.0, 
Vmax is the volume contained within 0.005 < z < 0.033). 
We assume that a LF does not evolve within the surveyed 
redshift range. 

We can test effectiveness of this volume correction by 
simulating galaxy distribution with a given LF. For sim- 
plicity, we adopt a flat LF such that galaxies are uniformly 



distributed within a magnitude range of —21 < M < —15. 
We examine a set of galaxy distribution here. 

In the top panels of Figure IAT1 we show a model of uni- 
form galaxy distribution in a given volume, and the middle 
of the top panels shows a volume corrected LF. As seen in 
the panel, a LF is reconstructed well. In the middle panels, 
we set non-uniform galaxy distribution in a given volume. If 
galaxies are uniformly distributed in a given volume, we ex- 
pect that the number of galaxies in a given radial shell scales 
as Ng a i oc r2 (this is the case for the top panels). But, here 
we adopt the distribution such that N ga i is independent of 
r. In this case, galaxies are not distributed uniformely in a 
volume, and a volume corrected LF does not reconstruct the 
parent LF. This is expected - the volume correction assumes 
the uniform galaxy distribution over a survey volume. The 
bottom panels show another example, where radial galaxy 
distribution has 'gaps'. The volume correction again fails to 
reproduce the parent LF. 

Therefore, we take into account variations in radial 
galaxy distribution such as large-scale structures when we 
perform the volume correction. In our sample, galaxies hav- 
ing My < —19.5 can be observed regardless of redshift (i.e., 
volume- limited) . We use these galaxies to correct for the 
radial variation. A number of galaxies with Mv < —19.5 
should scale as oc r2 if they are distributed uniformly in the 
Universe. We define radial distribution of galaxies having 
Mv < —19.5 divided by r2 as f(r), and radially averaged 
f(r) as fave- A statistical weight becomes: 

^i,corr — ( Vsurvey /Vmax ) X [/ave//(^)] - C^"^J 

Note that, in the first term, we do not include the volume 
where no galaxies are observed. For example, we do not cal- 
culate the volume contained in the 'gaps' in the bottom 
panels in Fig. IA1I because no statistical correction can be 
made where there are no galaxies. In the latter term, we 
make radial bins to compute f(r) and /„„. We confirmed 
that our conclusions are not strongly affected by uncertain- 
ties arising from the binning. The right panels of Fig. IA1I 
shows volume corrected LFs using this equation. The par- 
ent LFs is reconstructed well. Note that Eq. IA2I does not 
reproduce an absolute number of galaxies. 

As a final check of this correction, we use real galaxy 
distribution in our sample and artificially assign their mag- 
nitudes based on a flat LF. We confirm that the parent LF is 
recovered well in all the field, group, and cluster samples by 
using Eq. IA2I It turns out that the correction for the non- 
uniformity is particularly important for group and cluster 
galaxies whose distribution is strongly related to large-scale 
structures. Actually, cluster galaxies have 'gaps' in their 
radial distribution. To sum up, a flux-limited sample can 
be used to discuss galaxy properties by applying statistical 
weights of Eq. EH 



APPENDIX B: STATISTICAL 
CONTAMINATION SUBTRACTION 

This appendix describes the statistical contamination sub- 
traction procedure in detail. See Figure iBTl for help. We use 
the control field sample defined in J|]to subtract the con- 
tamination. For simplicity, in this appendix, we refer to the 
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Figure Al. The volume corrected LFs of three types of galaxy 
distribution. ^From top to bottom, the panels show uniform dis- 
tribution (N ga i oc r2), r-independent distribution [N„ a i is con- 
stant at a given r), and uniform distribution with redshift gaps. 
The left panels show galaxy distribution as a function of a dis- 
tance from us (unit is arbitrary). The middle panels present recon- 
structed LFs based on the Wi weighting. The right panels show 



reconstructed LFs based on the Wi 
details. 



weighting. See text for 
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Pick up a galaxy in the target sample 

Remember its magnitude and colour and move to the field CMD 

Pin down a point where the target galaxy should lie. 

Draw a circle (ellipse) around it 

Count field galaxies that fall in the circle. 

Repeat the above proceudre for all the target galaxies 

Then calculate field probabilities for the target galaxies defined by the ratio of 
the sum of statistical weights of field galaxies in the circle of the target galaxy 
to the total number of field galaxies 
Randomly pick up a target galaxy 

Throw a dice and determine if this galaxy is a contaminant galaxy or not according to 
its field probability. 
.Back to 8 and repeat 



Figure Bl. The contamination subtraction procedure. 



control field sample as field sample. The field sample ex- 
amined in the main text should be considered as a 'target' 
sample here. 

Let us pick up a galaxy in a target CMD. We move 
on to a CMD of field galaxies, and pin down a point where 
our target galaxy should lie (i.e., magnitude and colour of 
the galaxy we picked up). Then we draw a circle (strictly 
speaking, ellipse) around the point on the field CMD, and 
count field galaxies that fall in the circle. The above proce- 
dure is performed for all target galaxies, and then we know 
how many times a field galaxy is counted by the target 
galaxies. We define a statistical weight for a field galaxy 
as an inverse of its count. For example, if a field galaxy 
is counted three times (i.e., three target galaxies share the 
field galaxy), the statistical weight of the field galaxy is 1/3. 
Then we define a probability for a target galaxy to be a 
field galaxy as the ratio of the sum of statistical weights of 
field galaxies in the circle of the target galaxy to the to- 
tal number of field galaxies. We randomly pick up a target 
galaxy, throw a dice, and determine if the target galaxy is 
a field galaxy or not according to its field probability. We 
repeat this procedure until we subtract sufficient number of 
galaxies. In short, this procedure can be c onsidered as a gen- 
eralized form of the grid-ba sed method lIKodama fc Bowen 
l2001bl: IPimbblet et al1l2002f) . It should be noted that since 
we make no grids, we are free from an uncertainty on how we 
make the grids. Also note that we are free from the 'negative 
galaxy' problem which the grid-based method faces. 

There are two parameters in our method, namely, the 
radius and the shape of the aperture within which we count 
field galaxies. The aperture we actually use is shown in Fig- 
urc lBll It should be noted that our results do not depend 
on the size and shape of the aperture in reasonable ranges. 

How many galaxies should we subtract? When we have 



galaxies that occupy surface area of A arcmin2, we expect 
the contamination to be Eficid.4 galaxies, where Efi c id is an 
average surface density of field galaxies. Since we define envi- 
ronment by galaxy density, it is not very straightforward to 
estimate a surface area occupied by density-selected galax- 
ies. We infer a surface area by a simple approximation. Be- 
cause an inverse of a galaxy density means an average surface 
area occupied by a galaxy (i.e., unit of arcmin 2 galaxies^ 1 ), 
a sum of an inverse of density gives a surface area: 

N 1 

-4-= = E?--^, (Bl) 

l= 1 ' 

where Ei oca ( is local density and i runs through 1st to AT-th 
galaxy we select. We show in Appendix |C| that this esti- 
mate is statistically robust. A Poisson error is then added 
in Efieid^E in each Monte-Carlo realization of the contami- 
nation subtraction. 

We check if our statistical contamination subtraction 
technique properly subtracts field galaxies. We examine if 
subtracted target galaxies reproduce the CMD of field galax- 
ies. Figure lB"2l plots distribution of subtracted target galax- 
ies on a CMD projected onto 1 dimensional space for clarity. 
It can be seen that the distribution of subtracted galaxies 
closely follows that of field galaxies. The very small differ- 
ence between them is due to the fact that we apply a smooth- 
ing of the galaxy distribution on the field CMD when we es- 
timate field probabilities (the size and shape of the aperture 
determines the smoothing scale). Despite the small differ- 
ence we see a fairly good correspondence, and we therefore 
conclude that our subtraction method properly removes the 
field contamination. 
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Figure B2. Field galaxies in RXJ0153 are divided into 6 mag- 
nitude bins (6 panels from top to bottom), and fractions of field 
galaxies as a function of R — z' colour are plotted. The open 
squares show distribution of field galaxies. The solid circles show 
distribution of 'subtracted' target galaxies. 
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Figure CI. As / Aaper plotted against the number ratio of clus- 
ter galaxies to field galaxies in the aperture. The points mean 
the median of the distribution, and the error bars represent the 
quartiles of the distribution. 



APPENDIX D: GLOBAL DENSITY 



APPENDIX C: CONVERSION FROM DENSITY 
TO SURFACE AREA 

In Appendix [B] we have converted from local density to a 
surface area that is occupied by density-selected galaxies. 
We examine validity of this approximation in this section. 
We construct a simple toy model of galaxy cluster and in- 
vestigate differences between the real (metric) area and the 
area estimated from density. 

Our toy model consists of uniform field galaxies and 
cluster gala xies whose distribution follows the King model 
We randomly generate galaxies and estimate 
density for each galaxy by the nearest-neighbor method 
used in this paper. Galaxies within a fixed aperture around 
the cluster is selected. Then we compare the area of the 
aperture (Aaper) with the area estimated from densities of 
galaxies within the aperture {As , see eq. IB10 . We calculate 
As /Aaper in each realization and repeat 1000 times. There 
are two parameters in the model: the radius of the aperture, 
and the fraction of cluster galaxies to field galaxies in the 
aperture. 

As a fiducial value, the aperture radius is set to 3r c 
(r c means the core radius of the King model). Figure IHT1 
plots As /Aaper against the number ratio of cluster galax- 
ies to field galaxies in the aperture. The result is encourag- 
ing - As /Aaper is almost unity. In the density calculation, 
the distance to n-th nearest galaxy is used (n = 5 and 10 
is adopted in the main text). We run the simulation using 
n = 5, 10, 20, and find that there is a weak trend that the 
scatter in As /Aaper decreases with increasing n, but the 
median is always close to unity. Our choice of the aperture 
radius seems to have a rather strong effect. The scatter in 
As/ 'Aaper increases as we adopt smaller aperture (e.g., 2r c ). 
On the other hand, if we adopt larger aperture, say 5r c , the 
scatter reduces to almost half. However, in any case, the 
median As/ Aaper is always close to unity. Therefore, we 
conclude that As is a good measure in a statistical sense. 



Our aim in this appendix is to address how effective and 
quantitative our group/cluster separations are. For simplic- 
ity, groups and clusters are referred to as systems in this 
Appendix. 

We run the frie nds-of-friends algorithm (FOFA; 
iHuchra fc Gelled Il982h in our SDSS sample to find 
galaxy groups and clusters. The FOFA is a famous 
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FOFA is to find a set of galaxies connected within cer- 
tain linking lengths from each other. There are two linking 
lengths, namely, an angular separation (Do) and a line-of- 
sight velocity difference (Vb). Since galaxy systems are elon- 
gated along the line-of-sight in the redshift space ('finger of 
god' effect), we need to handle Do and Vb separately. 

We run the FOFA in a volume-limited sample defined 
by 0.010 < z < 0.065 and M g < -19.3 so that we do not 
have to apply any scaling in the linking lengths. The FOFA 
parameters are set to Do=700kpc and Vb = 400 km s -1 . 
These values are sli ghtly different from those adopted in 
iTanaka et alJ J2004l> . Our conclusions are not, however, 
strongly dependent on the choice of the parameters. There is 
one more parameter, N m i n , which determines the minimum 
size of galaxy systems. If a system has members less than 
Nmin, the system is not co nsidered here. We set N m i n — 5 
as m ITanaka et aP J2004T) since a signific ant fraction of 
N < 5 groups are expected to be spurious <Frederidll99i 
IRamella. P isani. & Gelledll997t iRamella et al.ll2002l) . 

In Figure ITjTl we present a correlation between the num- 
ber of member galaxies (N ga i) of the FOFA systems and 
global density of the member galaxies. The correlation is 
encouragingly good - there is a positive correlation between 
the two quantities, especially at N ga i > 20. N ga i < 20 sys- 
tems show only a weak correlation. This is probably because 
that a typical extent of such poor systems is smaller than 
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Figure Dl. Global density plotted against N ga i determined 
by the FOFA. The dot and the error bars show the median and 
quartiles of global density distribution in each system. 

2Mpc and a number fluctuation of field galaxies contribute 
to global density. As for N ga i > 20, cluster members dom- 
inate the 2Mpc aperture and hence we have a good corre- 
lation there. A few systems that have small N ga i with high 
global density can be seen. They are either parts of rich 
systems accidentally splitted by the FOFA or real isolated 
systems close to nearby rich systems. They comprise only a 
small fraction and the effect of such contamination is negligi- 
ble. In summary, global density is a powerful tool to separate 
groups from clusters quantitatively. 

We measure velocity dispe rsions of detected systems 
using the biweight estimator feeers. Flvnn. fc Gebhardl 
1990). By examining the relationship between N ga i and 
a velocity dispersion of a system (a), we find that the 
global density threshold used in the main text for the 
SDSS, Egiobai = 2.5, roughly corresponds to systems having 
a = 300 - 400kms _1 . Therefore, our groups typically have 
a < 300-400 kms -1 and clusters have a > 300-400 kms" 1 . 



